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ABSTRACT 

Many models of gamma-ray bursts (GRBs) involve accretion onto a compact object, 
usually a black hole, at a mass accretion rate of order a fraction of a solar mass per 
second. If the accretion disk is larger than a few tens or hundreds of Schwarzschild 
radii, the accretion will proceed via a convection-dominated accretion flow (CDAF) in 
which most of the matter escapes to infinity rather than falling onto the black hole. 
Models involving the mergers of black hole white dwarf binaries and black hole helium 
star binaries fall in this category. These models are unlikely to produce GRBs since 
very little mass reaches the black hole. If the accretion disk is smaller, then accretion 
will proceed via neutrino cooling in a neutrino-dominated accretion disk (NDAF) and 
most of the mass will reach the center. Models involving the mergers of double neutron 
star binaries and black hole neutron star binaries fall in this category and are capable 
of producing bright GRBs. If the viscosity parameter a in the NDAF has a standard 
value ~ 0.1, these mergers can explain short GRBs with durations under a second, but 
they are unlikely to produce long GRBs with durations of tens or hundred of seconds. 
If the accretion disk is fed by fallback of material after a supernova explosion, as in 
the collapsar model, then the time scale of the burst is determined by fallback, not 
accretion. Such a model can produce long GRBs. Fallback models again require that 
the accretion should proceed via an NDAF rather than a CDAF in order for a significant 
amount of mass to reach the black hole. This condition imposes an upper limit on the 
radius of injection of the gas. 



1. Introduction 

The fireball model (see Piran 1999, 2000, for reviews) provides a good understanding of 
conditions within the 7-ray-emitting and afterglow-emitting regions of gamma-ray bursts (GRBs). 
According to this model, GRBs are produced when relativistic ejecta from a "central engine" 
are slowed down by interactions, either with an external medium (the external shock model) or 
among different layers within the ejecta themselves (the internal shock model). In the interactions 
the kinetic energy in the ejecta is converted to relativistic electrons which produce the observed 
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radiation. Among the many successes of the model we note the observational confirmation of 
relativistic motion in the afterglow (Frail et al. 1997; Katz & Piran 1997). 

Despite the successes of the fireball model, the nature of the central engine remains a mystery. 
The problem is that the central engine is hidden from view; no radiation (apart from gravitational 
radiation and neutrinos that may possibly be detected in the distant future) reaches the observer 
directly from the engine. For a number of so-called "long bursts" accurate positions have been 
determined through observations of their afterglows. Based on this, there is circumstantial evidence 
that these bursts are associated with star-forming regions (e.g. Bloom et al. 2000). There is no 
information at present on the other class of GRBs, the so-called "short bursts." 

Although we lack direct evidence on the nature of the central engine, it is nevertheless widely 
accepted that GRBs are the result of cataclysmic events involving either neutron stars or stellar- 
mass black holes. The arguments in support of this hypothesis arc straightforward, (i) Since bursts 
radiate the bulk of their energy in the 7-ray band, it seems likely that a relativistic object is behind 
their production, (ii) The energy budget (~ 10^^ erg) is comparable to the kinetic energy of ejecta 
in a supernova explosion, (iii) Most long bursts are highly variable in gamma-rays, and so are 
many short bursts (Nakar k, Piran, 2001a). In particular, the ratio of the total duration of the 
burst to the variability time scale is large, from which one concludes that the gamma-rays must 
be produced in internal shocks (Sari & Piran, 1997). A key feature of internal shocks is that the 
observed gamma-ray variability reflects the variability in the activity of the central engine (Sari &: 
Piran, 1997). Since variability time scales as short as a millisecond are observed, the engine must 
contain a compact object of no more than a few solar masses (otherwise the light-crossing time 
would exceed the variability time). 

An interesting clue to the nature of GRBs is provided by the durations of bursts. While 
the fastest variability time scale is under a millisecond, burst durations are usually very much 
longer. Long bursts have durations ranging from 10-1000 seconds, and even short bursts have a 
median duration of about 0.3 seconds. Clearly, whatever is the physical mechanism behind GRB 
production, it acts on a much longer time scale than the fastest dynamical time of the central 
engine. 

Narayan, Paczyhski & Piran (1992) suggested that the central engine in GRBs involves the 
accretion of matter onto a compact star, and that the energy in the burst is provided by the 
gravitational energy released by the accreting gas. In such a model, the duration of the burst is 
set by the viscous time scale of the accreting gas. In most accretion flows, the viscous time is 
significantly longer than the dynamical time, and so the accretion model naturally explains the 
large difference between the durations of bursts and their fastest variability time. 

The formation of an accretion disk is a natural outcome of most popular models of GRBs, 
e.g. the mergers of double neutron star binaries (Eichler et al., 1989; Narayan, Paczyhski & Piran 
1992), neutron star black hole binaries (Paczyhski, 1991; Narayan Paczyhski & Piran 1992), black 
hole white dwarf binaries (Fryer et al., 1999), black hole helium star binaries (Fryer & Woosley, 
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1998), and models based on "failed supernovae" or "collapsars" (Woosley, 1993; Paczynski, 1998; 
MacFadyen & Woosley, 1999). An important exception is Usov's (1992) model in which the GRB 
energy is provided by the magnetic and rotational energy of a newly formed rapidly rotating neutron 
star. 

In this paper, we consider a generic accretion model of a GRB in which a certain amount 
of mass, mdMQ, goes into orbit around a relativistic star of mass Stji^Mq. We assume that the 
orbiting mass is initially inserted into a torus at a radius routes , where Rs is the Schwarzschild 
radius of the central star: Rs = 2GM/(? = 8.85 x 10^ m3 cm. Starting from the initial toroidal 
configuration, the mass spreads out by viscosity and becomes an accretion flow extending from 
r = 1 (the horizon of the central black hole) to r ~ rout- (There might also be an outflow at 
radii > rout as we discuss below.) We work out in this paper the time scale, tacc, of the accretion 
flow, the average mass accretion rate, M = mMQS~^, onto the central star, the amount of mass, 
fna.cc = iaccfn, accreted by the star, and the accretion efficiency, ^ = macc/^n^d- 

These parameters are constrained by observations. In binary merger models, the durations of 
bursts should be comparable to the accretion time iacc- Durations are of order a second or less for 
the class of "short duration GRBs" and in the range 10-1000 seconds for the class of "long duration 
GRBs." The energy in long duration bursts is estimated to be ~ 5 x 10^° erg (Panaitescu &; Kumar, 
2001; Frail et al. 2001). With a reasonable efficiency of converting accretion energy to relativistic 
flow (< 0.01) this corresponds to rriacc ~ O.IMq and to a peak accretion rate M > IO^^A/qS^"*^. The 
values of these parameters for short bursts are less certain, as the distance scale to this population 
is uncertain. However, it is unlikely that they are smaller by more than an order of magnitude. 

The mass accretion rate tends to be extremely high for a typical GRB model — it is of order 
a fraction of a solar mass per second. At such accretion rates, the optical depth of the accreting 
gas is enormous and radiation is trapped inside the gas. In the normal course, the accretion would 
proceed via a radiatively inefficient flow, such as an advection-dominated accretion flow or the 
related convection-dominated accretion flow. We consider such flows in §2. If the accreting gas 
has a sufficiently high temperature and density, however, it can cool via neutrino emission, leading 
to a neutrino-dominated accretion flow, as discussed by Popham, Woosley & Fryer (1999) (see 
also Ruffert k, Janka, 1999, for a numerical simulation of a binary merger which included neutrino 
losses). We discuss this case in §3. In §4 we go beyond the analytical results of §§2,3 and present 
numerical results, delineating the regions of (rout) ?^d) space where different types of accretion 
occur. We discuss in §5 the implications of the results for various models of GRBs. 



2. Radiatively Inefficient Accretion: CDAF 

At the mass accretion rates of interest to us, the optical depth of the gas is extremely high, 
and the radiation is very effectively trapped (we demonstrate this below). The accretion flow 
then corresponds to an advection-dominated accretion flow, or ADAF (Narayan &; Yi 1994, 1995; 
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Abramowicz et al. 1995; Chen et al. 1995; see Narayan, Mahadevan & Quataert 1998 and Kato, 
Fukue & Mineshige 1998 for reviews). Radiation-trapped ADAFs were initially discussed by Katz 
(1977) and Begelman (1978), and later analyzed via height-integrated "slim disk equations" by 
Abramowicz et al. (1988). 

Igumcnshchcv and co-workers have carried out hydrodynamic simulations of ADAFs and have 
discovered several interesting properties of these flows (Igumenshchev, Chen & Abramowicz 1996; 
Igumenshchev & Abramowicz 1999, 2000). They find that, when the dimensionless viscosity pa- 
rameter a is large, say a. > 0.3, the accretion flow has a strong bipolar outflow, as anticipated in 
some previous papers (Narayan h Yi 1994, 1995; Blandford h Begelman 1999). On the other hand, 
when the viscosity is relatively weak, say a < 0.1, the flow has well-developed convection (which 
was again anticipated in earlier work, cf. Narayan & Yi 1994, 1995). A convective ADAF has quite 
unusual properties (Stone, Pringle & Begelman 1999; Narayan, Igumenshchev & Abramowicz 2000; 
Quataert h Gruzinov 2000; Igumenshchev, Abramowicz h Narayan 2000), which arise because con- 
vection moves angular momentum inward rather than outward (Narayan et al. 2000; Quataert &; 
Gruzinov 2000). Convective ADAFs have been given the name of convection-dominated accretion 
flows, or CDAFs. 

As in other accretion flows, "viscosity" in radiation-trapped flows is believed to arise from 
magnetic stresses resulting from the Balbus-Hawley instability (Balbus & Hawley, 1991). Numerical 
simulations of shearing MHD flows give values of a in the range a. < 0.1. We therefore assume that 
radiation-trapped accretion flows in GRBs also have a in this range (no 3D MHD simulations of 
these flows have been done so far, but there has been some 2D work by Stone h Pringle 2001). 
Given the relatively low value of a, we expect the accretion flow to take the form of a CDAF. 

Ball, Narayan h Quataert (2001) have given approximate scalings for various fluid variables 
in a CDAF. The density and velocity scale as 

p{r) « Poutf— ) ' =2.97x lO^^mgWo^uA''/' gcm-3, 

\ ''out / 

v{t) ^ cr-=^/2 = 3 X 10lV-=^/^ (1) 

where, in the first equation, we have used 4Trr^^^R^{H/R)pont = n^d^Q to relate pout in the CDAF 
to the mass in the accretion flow: 

/9out = 2.97 X lOi^mg^mdr-^ g ^^^-3 (2) 

Note the use of in rather than 47r/3 in the equation for the mass. This has been done to obtain a 
better match with the results for an NDAF (see §3). 

We assume that the accreting gas consists of photo-disintegrated nuclei with roughly equal 
numbers of neutrons and protons (and electrons). The optical depth through the flow is 

Tout ~ ^^cTTVoutRs = 5.21 X 10^^ m^'^ mar (3) 
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The diffusion time for radiation to leak out of the flow is then 

idiff ~ Tout^^^^^ = 1.54 X lOi^mJ^dr-^t s. (4) 

This time is very much longer than the accretion time for the parameter ranges of interest to us, 
namely ~ 1, ~ 0.1 — 1, rout ~ 10 — 10^. Thus we expect radiation to be very effectively 
trapped within the accretion flow. 

Ball et al. (2001) also estimate the isothermal sound speed Cg and the scale height H of 
the accretion flow. We use slightly different coefficients here which are more appropriate for a 
radiation-pressure dominated 7 = 4/3 gas (as opposed to the 7 = 5/3 gas that Ball et al. 2001, 
considered) : 

cl ^ 0.3cV-^ = 2.70 X lO^^r-^ cm^ s'^ 

These scalings are approximate, but they are likely to be accurate enough for the purposes of this 
paper. The mass accretion rate in the CDAF is estimated to be M = Airr'^ Rg{H / R) p{r)v{r) , which 
gives, using the above relation for H/R, 

m = 3.39 X 10"^ m^^mdr'^/'^. (6) 

(Recall that rh is defined as the mass accretion rate in units of solar masses per second.) 

The accretion time scale is not simply equal to md/rh. The reason is that much of the mass in 
a CDAF actually flows out of the system rather than into the central black hole (Stone et al. 1999; 
Narayan et al. 2000; Quataert &; Gruzinov 2000Igumenshchev &: Abramowicz 2000). We therefore 
proceed as follows. 

The random velocities of convective blobs in a CDAF are typically less than the local Keplerian 
velocity by a factor which is proportional to the viscosity parameter a (see Narayan et al. 2000). 
We thus write approximately 

^^turb ~ C(Vk = 2.12 X 10^a_ir-^/2 cms'^ (7) 
The residence time of a convective blob at radius r is then 

ires ~ = 4.17 X IQ-^aZlmsr^/^ s. (8) 
We make the reasonable assumption that the accretion time is of the same order as ires at r = rout^ 

iacc ~ *res(r-out) = 4.17 X IQ-^aZlmsrli^^ s. (9) 

We then find that the amount of mass accreted by the black hole is mace = ^acc'^) but since this 
mass cannot exceed the total available mass irid we write: 



ma 



rrid, rout < U.la_lmd, 

14.lQ;lJm<irout, ^out > UAaZlmd- (10) 
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Note that, when rout is large, the accreted mass is much less than m^. The reason is that the 
bulk of the mass is ejected from the system, flowing out at r ~ ^out- The energy for the ejection is 
provided by convective energy flux from the interior of the flow. 

The coefficient 14 in equation (10) is somewhat uncertain since wc do not know the exact 
relation between tacc and t^cs] there could well be a numerical factor other than unity relating the 
two. A different (and more more detailed) derivation of equation (10) is given in the Appendix, 
where we again find that the coefficient is uncertain. In the rest of the paper we use equation (10) 
as written, but we should keep in mind that it could be in error by a factor of a few. 

Let us now calculate the temperature of the CDAF. The pressure in the accreting gas is given 

by 

p = pcl = 8.02 X lO^'^m^^mdrJ/'^r-^/^ ergcm^^ (11) 

The pressure has three contributions: radiation pressure, gas pressure, and degeneracy pressure 
(Popham et al. 1999): 

11 4 pkT 2iThc f 3 \^^^ { p\^^^ 
^12 rup 3 \87rmp/ \Me/ ^ ^ 

The quantity a is the radiation constant, and the factor 11/12 includes the contribution of relativis- 
tic electron-positron pairs (asssuming that the temperature is sufficiently above the pair threshold 
limit). The gas pressure term includes the contributions from non-relativistic particles. We assume 
that we have an equal mix of protons and neutrons (i.e. we assume that all complex nuclei have 
been photo-disintegrated). In the degeneracy pressure term we use an electron molecular weight 
jjLe = 2, assuming an equal mix of protons and neutrons; the ratio of neutrons to protons in the 
interior of a NS is about 0.2. The contribution of e~-e^ pairs has been ignored in the degeneracy 
pressure term. For the calculations presented in §4, we solve for T by using the full expression for 
p given in equation (12). But here, in order to obtain simple analytic estimates, we simplify the 
equation and assume that radiation pressure dominates. We then find that 

T = 1.84 X lO^^m- '/'myV-f r-3/« K. (13) 

This estimate of T is valid only if gas pressure is smaller than radiation pressure. The gas pressure 
is equal to 

p^^ = Pj^= 4.52 X lO^^m- ''/"m^V^-f /^-^/« ergcm-=^. (14) 

Vflp 

For this to be smaller than radiation pressure we require 



r 



< l.hm^J^rrq'^'^. (15) 



^'out 

For our fiducial black hole of mass SM©, gas pressure dominates only if the disk is quite massive, 
rrid ^ 5, and if r is close to rout- (The detailed numerical calculations in §4 show that gas pressure 
dominates for ma ^ 1 rather than 5, but even for such values of m^, if r is much less than rout; 
radiation pressure takes over.) 
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Let us next calculate the rate of cooling of the accreting gas. We showed earlier that cooling 
via radiative diffusion is negligible. However, at the temperatures found in these flows, cooling via 
neutrino losses may be important (cf. Popham et al. 1999). The cooling rate per unit volume due 
to neutrinos, , takes the form 

Q~ = QuD + ^Jiv ~ 5 X lO^^Tfi + 9.0 x lO^V^n erg cm'^ s-\ (16) 

where the first term on the right describes cooling via pair annihilation (the so-called URCA 
process) and the second term describes cooling via pair capture on nuclei (estimated for Xnuc = 1 
as appropriate for our fully photo-disintegrated nuclear gas, cf. Popham et al. 1999). For the 
range of parameters of interest to us, invariably dominates over q~p. We therefore neglect q~p 
hereafter. 

We then obtain the following estimate for the cooling rate in the CDAF: 

q^j, = 1.04 X lO^^m-^'/'mf r-f ergcm-^g-i. (17) 

Since the energy density of the radiation-dominated gas is equal to 3prad = (11/4)0^^, we estimate 
the cooling time of the gas at radius r to be 

tcooi = ^ = 2.31 X 10--m^/^m-^/V2V^/^ s. (18) 

leN 

One of the important requirements for the accretion flow to behave like a CDAF is that it should 
be radiativcly incfRcicnt. We, therefore, require the cooling time tcool at each radius to be longer 
than the residence time tres of a convecting blob at that radius. The ratio of cooling to residence 
time is 

^ = 5.54 X 10-«a_impm;'/VSV-V4. (19) 

Since the ratio decreases with increasing r, we set r = rout- Then we see that a CDAF is possible 
only if the following condition is satisfied: 

rout > llSali/^mJ^mJ/^. (20) 

If the mass in the accretion flow is initially at a radius greater than the above limit, then the flow 
will become a CDAF at r = rout; and indeed for all r < rout- The scaling relations derived in this 
section would then be valid. If, however, the initial radius of the gas is below the above limit, then 
neutrino-cooling will prevail, and we will have a cooling-dominated accretion flow at all radii from 

rout down to the black hole. 

For completeness, we consider also the case when gas pressure dominates (which requires a 
large value of nid, and r close to rout) as shown above). The temperature is then obtained by setting 
P = Pgas- This gives 

T = ^ = 3.27 X lO^V-i K. (21) 
k 
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The analysis presented so far assumes that we are given the initial mass in the disk and the 
initial radius of the gas rout- Another possibility is that we have steady injection of mass at a rate 
minjM0S~^ at a circularization radius rout- This is the case, for instance, in the collapsar model, 
where the material is supplied by fallback from a supernova explosion. We may assume that the 
accretion flow achieves an approximate steady state in which the rate of injection of mass equals 
the rate of mass loss from the flow (both inflow and outflow): minj = ?TZd/iacc- Therefore, we obtain 

ma = 4.17 x W~'^aZ\m3mnirlil (22) 

This expression, which gives the mapping between rrid and rhinj, may be substituted into the 
various relations derived in this section to obtain the corresponding results for the case of steady 
mass injection. 



3. Radiatively Efficient Accretion: NDAF 

We saw in the previous section that when the mass is initially injected at an outer radius 
smaller than the limit given in equation (20), cooling via neutrino emission becomes significant 
and the flow is no longer radiatively inefficient. We then have a regime of accretion in which the 
viscous energy dissipation is balanced by neutrino cooling. Popham et al. (1999) named this a 
neutrino-dominated accretion flow (NDAF) and worked out its properties. We extend their results 
in this section. 

Because the gas cools efficiently, an NDAF behaves like a thin accretion disk and we are entitled 
to use the basic theory of thin disks (Shakura Sz Sunyaev 1973; Frank, King & Raine 1992). As 
before, the pressure of the gas has three contributions and is described by equation (12). The 
isothermal sound speed Cg and the vertical scale height H are given hy = p/ p and H = Cs/^k, 
where = {GM/R^)^'^ = 2 AO X lO^mg s ^ is the Keplerian angular velocity. We write 

the coefficient of kinematic viscosity in the usual form as 

^ = a-f, (23) 

where a is a dimensionless parameter which we expect to have a value 0.1. We denote, therefore, 
a-i = a/0.1, and write down our scalings in terms of a-i. 

A thin accretion disk satisfies two equations. Angular momentum balance gives 

R ^ -1/2] 



1 



R* 



GTTupH, (24) 



where T, = 2pH is the surface density and R^ is the radius of the inner edge of the disk (3-^5 for a 
Schwarzschild black hole) . The approximation on the right is valid for R^ R^,. The radial velocity 
of the gas is given by v = Su/IR and so the accretion time is 

^ _ -Rout _ 2-Rout ^2^-^ 



-9- 



If we write the disk mass as rrid = ts^M, then by combining the previous two equations we find 
that 

md = 4.'KRl^^(-\ Pont- (26) 
^/ ( 



out 

Note that we used this relation, with the same coefficient Att, also for a CDAF. 

The condition of energy balance (viscous heating equals radiative losses) gives: 

3GMM , _ 

q H = {q^j,+q,^)H. (27) 



This relation closes our set of equations and allows us to solve for p, T and other quantities. We 
show numerical solutions of the full equations in §4. 

In the rest of this section we derive analytical scalings by making some approximations. First, 
the numerical calculations show that the g^^y term in the cooling law dominates over q^^j for all 
parameters of interest. We will therefore assume this. The calculations also show that for most 
parameters, either gas pressure or degeneracy pressure dominates. (There is a small region of 
parameter space where radiation pressure dominates; we ignore this region in the analytical work.) 

Let us first assume that gas pressure dominates. Expressing all results in terms of the scaled 
temperature, Tn = T/lO^^if, and substituting in the angular momentum equation, we obtain a 
relation between p and Tn 

p = 2.56 X 10^^aZ\m^'^mr-^T~^^^'^ gcm"^. (28) 

Substituting this in the energy equation, we can solve for Tn, and thereby obtain the various other 
quantities (Popham et al. have derived similar relations): 

W = 2.76 X lO-^alf mf r^S, 
rh = 36.2a^_/^ 771^*^^^ nidr^^/^, (29) 

o oo inl5 -1/10 -29/10 -67/20 -3 

Pont = 2.28 X 10'°a_/ ' rudv^^^' gem 

In a cooling-dominated thin disk, very little mass is expected to be lost to outflows, so we expect 
nearly all the mass in the disk to be accreted by the star, i.e. 

rrie^c ~ rud. (30) 

The above results are valid provided gas pressure dominates. By comparing gas pressure to 
degeneracy pressure we can determine the condition for this to be true. We find that we require 

Tont > 26.2a-_rm-''/''mf''. (31) 
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When the outer radius is below this hmit, degeneracy pressure takes over from gas pressure and 
we obtain a different set of analytical scalings: 

. orr -13/7 9/7 -3/2 

Pout = 7.32 X lO^^m- ^'/^/Vut 

Tn = QmQa'ltrn-''''^mf\:T\ (32) 
W = 2.82 X lO-^alJmf/^m-'/V^/t s. 

In our analysis of both the gas-pressure-dominated and degeneracy-pressure-dominated regimes, 
we assumed that the accreting gas is optically thin to its own neutrino emission. This assumption 
breaks down at sufficiently small radii. We may estimate the neutrino optical depth as 

r. = = 5.20 X l^'aZrrnt''''n^T'r-^T■ (33) 

The radius at which the optical depth goes to unity is 

r.=i = 2Q.ha-T'm-''"'''mf\ (34) 

This radius lies within the degeneracy-pressure-dominated zone. Inside this radius, we need to 
consider neutrino transport in more detail. This is beyond the scope of the paper. 

It is interesting to examine whether the NDAF solution is stable. Following Piran (1978) we 
use the general condition for thermal stability: 

/dlnQ+\ ( d\nQ- \ 

where are the integrated (over the height of the disk) heating (-I-) and cooling (— ) rates. For 
an NDAF with viscosity described by the a-prescription, the heating rate goes as Q"*" oc i?^, and 

the vertically integrated cooling rate (for pair capture on nuclei) goes as Q- oc T.T^. For an NDAF 
in which gas pressure dominates, T oc H'^ and oc T,H^^. The criterion (35) is satisfied easily 
and we conclude that such an NDAF is thermally stable. 

If radiation pressure dominates, T oc p oc Si^ and Q- oc Y.^/'^H^I'^. The condition (35) is not 
satisfied and the NDAF is unstable. This resembles the situation in "conventional" accretion disks 
whose inner regions become thermally unstable when radiation pressure dominates (Frank et al. 
1992). 

If degeneracy pressure dominates, the temperature is independent of either or S. The 
temperature of the disk can then freely adjust such that Q~ balances any . The situation is 
clearly thermally stable. 



The condition for viscous stability is: 

dM 



ffi > 0. (36) 
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We have M oc S for the gas pressure case, M oc for the radiation pressure case, and M oc S^/^ 
for the degeneracy pressure case. All three cases are viscously stable. 

Thus, combining all these results on stability, we find that an optically thin NADF is unstable 
only if it is radiation pressure dominated. There is a very narrow region of parameter space near the 
boundary between NDAFs and CDAFs in (routj space where radiation pressure does dominate. 
This unstable region could conceivably play a role in determining the temporal behavior of some 
bursts. The stability properties of the optically thick regions of the NDAF remain to be worked 
out. 

The disks that we are considering are massive and we should also consider gravitational insta- 
bilities. The CDAF zone is always gravitationally stable in our models; the Toomre Q parameter is 
invariably much greater than unity. For gas pressure dominated NDAFs, we find that the Toomre 
Q parameter is given by 

Q = 8.5 a^rruTWr-^'^'r^J^t- (37) 

We see that Q decreases with increasing r, so that the flow is most unstable on the outside. However, 
even for r = Tom*, Q is sufficiently larger than unity (for all reasonable values of parameters) that 
we arc guaranteed stability. A similar result applies to radiation pressure dominated NDAFs. 

The case of degeneracy pressure dominated NDAFs is more complicated. For most of the 
parameter space {rrid, Tout) we find Q to be greater than unity, but it is only marginally so (Q ~ 3). 
When Tout < 10, we find that Q might become less than unity, signifying gravitational instability. 
But this is also the region in which the disk becomes optically thick to neutrinos and the analysis 
we have carried out breaks down. 

4. Numerical Results 

We discuss in this section numerical results which we have obtained by solving the full equa- 
tions. We assumed that we are given the initial mass of the accretion disk and the initial radius 
Tout- For each choice of rout and m^, we first assumed that the flow consists of a CDAF and used 
the relations described in §2 to calculate the properties of the accreting gas. In particular, we used 
equations (11) and (12), with p given by equation (1), to solve for the temperature T as a function 
of r. We checked the flow at all radii from r = rout down to r = 1 to make sure that the gas is 
radiatively inefficient at all radii. Specifically, we estimated the cooling time tcoob using the cooling 
formula given in equation (16), and checked that tcool is longer than the residence time of convective 
blobs ties at that radius. If tcool > ^rcs at all r, then wc identified the flow as a pure CDAF and 
estimated the accretion time tacc and the total mass accreted maco using the formulae given in §2. 

For a small region of parameter space near the transition between a pure CDAF and a pure 
NDAF, we found that the flow starts off as a CDAF at r = rout and switches to an NDAF at a 
smaller radius. We estimated tacc and rriacc for these cases by using a reasonable matching formula 
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at the transition radius. (We do not provide the details here, since this case is seen for only a small 
range of parameters). 

For Tout less than a critical value (whose value depends on rrid), we found that the flow cools 
too rapidly at r = rout to be a CDAF. In such cases, the flow becomes an NDAF at r = rout and 
remains an NDAF all the way down to r = 1. We solved the corresponding set of equations (see 
§3), and estimated tacc and mace appropriately. 

Figures 1 & 2 show contours of tacc and the accretion "efficiency" ^ = rrin^cc/^d, plotted in 
the space of the two principal parameters of the problem, rout and rrid- These models have the 
"canonical values" = a_i = 1. The boundary between the CDAF and NDAF zones is clearly 
seen in both figures, but especially in Fig. 1. The numerically determined location of the CDAF- 
NDAF boundary is fairly close to the analytical approximation given in equation (20) and shown 
as the right-most dotted line in Fig. 1. 

As equation (9) shows, the accretion time scale tacc for a CDAF is independent of rrid and 
depends only on rout- In contrast, tacc has a more complicated dependence on rout and nid for an 
NDAF (eqs. 29 and 32). These behaviors are seen clearly in Fig. 1. 

For a given initial mass m^, the accreted mass mace is equal to nid, independent of rout; in 
the case of an NDAF, and so ^ = 1 (see Fig. 2). However, if the flow becomes a CDAF (which 
happens for rout greater than a critical value), mace is signiflcantly less than nid, and the rest of the 
mass is ejected from the system; mace 

scales as rout i^ this case (cf eq 10). Therefore, ^ becomes 
significantly less than unity. Figure 2 illustrates this dependence. 

In some GRB models, such as those involving a BH-NS merger or the collapse of a very massive 
star, the mass of the black hole could be larger than 3Mq. We have therefore computed results 
for an M = 3OM0 black hole, i.e. ms = 10. Equations (20), (31) and (34) show that the various 
critical radii vary roughly inversely as (which is equivalent to saying that the physical radii 
R = rRs at which the corresponding transitions occur are independent of the black hole mass). 
We have confirmed this result in the detailed numerical calculations. As a result, for ra^ = 10, the 
NDAF zone shrinks by a factor of 10 in rout; and the various variants of the NDAF (degeneracy 
pressure dominated zone, optically thick zone) practically disappear. Thus, large black hole masses 
are not conducive to the formation of an interesting NDAF zone. 

In addition, we investigated models with a smaller value of the viscosity coefficient: a = 0.01, 
i.e. a-i = 0.1. In this case, the various critical radii become larger by roughly a factor of 2, as 
expected from equations (20), (31) and (34). More importantly, tacc increases by a factor of ~ 10 
(see equation 25). 

We also considered the case when the accretion is fed at a constant rate minj (rather than being 
initiated with an instantaneous addition of mass m^ as we have assumed so far). This case could 
be relevant for collapsar models. Figure 3 shows contours of accretion cfl&cicncy defined in this case 
as ^ = m/minj, plotted in the space of rout and rhinj for the "canonical values" ma = q_i = 1. The 
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overall behavior is consistent with the scaling relations derived in §§2,3. 

5. Discussion 

The starting point for this work is the fact that the accretion flow in a putative GRB central 

engine can have two very different forms: we expect the accretion to occur as a radiatively inefficient 
CDAF (convection-dominated accretion flow, cf Narayan et al. 2000, Quataert &: Gruzinov 2000) 
if the mass is introduced at a somewhat large outer radius, namely rout greater than the limit 
given in equation (20), while we expect the accretion to proceed via a radiatively efficient NDAF 
(neutrino-dominated accretion flow, cf. Popham et al. 1999) if rout is smaller than this limit. For 
a narrow zone in parameter space close to the CDAF/NDAF transition, it is possible for the flow 
to be a CDAF on the outside (r < rout) and to switch to an NDAF on the inside. But this is rare. 
By and large, for most choices of rout and m^^, the flow is either a CDAF at all radii or an NDAF 
at all radii. Since the two kinds of flow are very different from each other, there is a rather large 
difference in what an observer would see in the two cases. 

In most of the region of (routiiT^d) space where the flow is a CDAF, the mass accretion rate m 
and the amount of mass accreted mace are both very small. This is because in a CDAF, especially 
when rout is large, much more mass flows out of the system than into the black hole. If, as seems 
reasonable, a GRB engine requires a relatively large rriacc in order to produce a viable burst, then 
our work suggests that systems which form CDAFs are less likely to produce observable bursts. 

In contrast, an NDAF has substantially larger values of m and mace for a given m^. (In fact, 
we assume that mace = "md for an NDAF; there is probably some mass loss in a wind, but we expect 
it to be a small fraction of the total mass.) Therefore, a source with an NDAF is a much more 
plausible model of a GRB engine. By this argument, neutrino cooling (the key ingredient of an 
NDAF) is important for an efficient GRB. Note that we are not assuming anything about the actual 
mechanism of a burst. The fireball may be produced through neutrino-antineutrino annihilation 
(Eichler et al., 1989), or it could be the result of some other mechanism. Our suggestion is that 
whatever the mechanism may be, if it is based on accretion, then it probably requires a large mace 
to operate efficiently; such an mace can be achieved only through neutrino cooling of the accreting 
gas in an NDAF model. 

We note a possible caveat to the above conclusion. The energy in long duration GRBs is 
estimated to be about 5 x lO^'^ erg (Panaitescu &: Kumar, 2001; Frail et al. 2001). This energy 
could in principle be generated in a CDAF model with mace ~ O.IM©, provided the gravitational 
energy of the accreted matter is used very efficiently to launch a high Lorentz factor wind. We 
feel, however, that it is reasonable to hypothesize that the efficiency for converting the gravitational 
energy of maec to a GRB-producing relativistic wind is low, say not larger than 1%, and that bursts 
arise only from systems that form NDAFs and have large maec- 

Figures 1 and 2 show results corresponding to a simple model in which we assume that the 
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accretion disk is formed instantaneously, e.g. by the disruption of a companion star. We specify 

the initial state of the disk by giving its mass and initial radius rout- For such a model, we see 
that the region of parameter space where NDAFs form corresponds to short accretion times ^acc 
of order a few tenths of a second. This suggests that a binary-disruption-based accretion model is 
capable of producing short GRBs. It is, however, very hard to see how such an accretion system 
could produce a long burst. 

One way to make a long burst from an NDAF is to decrease the value of the viscosity parameter 
a (see eqs. 29 and 32). Could a be significantly lower, say 0.01? MHD simulations of thin accretion 
disks generally give a in the range 0.01 upwards, and it is widely agreed that the values obtained 
are lower limits since the simulations have limited spatial resolution. Empirical estimates of a 
in cataclysmic variables (obtained by comparing observations of dwarf nova outbursts with model 
predictions) give a ~ 0.1 in the high state (which is most relevant for our models). There is clear 
evidence that a is smaller, perhaps ~ 0.01, in the low state of CVs. This is probably the result 
of the cold gas becoming neutral and losing its coupling to the magnetic field (Gammie & Menou 

1998) , which is clearly not relevant for our ultra- hot plasma. We feel that a ~ 0.1 is a reasonable 
estimate for the viscosity parameter. 

Sakimoto &; Coroniti (1981) suggested that a may be lower in radiation-dominated gases 
because the magnetic pressure may be in equipartition with only the gas pressure rather than the 
total pressure. Their proposal requires that the gas and the radiation be able to slip past each 
other (cf. Blaes &; Socrates 2001). At the extraordinarily large radiation optical depths found in 
our models (both CDAFs and NDAFs), the gas and the radiation are extremely tightly coupled. 
It is therefore very unlikely that the value of a would be modified. 

The models we consider have large disk masses. The accreting gas may therefore develop 
gravitational instabilities (see the discussion near the end of §3) and lose angular momentum via 
gravitational waves (see Bonnell & Pringle 1995 and references therein). This might lead to an 
increase in the effective value of a. Another potential Tincertainty is that, for some choices of the 
parameters, electron degeneracy becomes important. It is not understood how the Balbus-Hawley 
instability (which is thought to produce the shear stresses behind a) behaves in such a gas. 

While the time scale of the accretion tacc depends fairly sensitively on a, the size of the NDAF 
zone is insensitive to a (cf eq. 20). Thus, Figs. 1 and 2 give fairly reliable limits on the radius inside 
which the mass needs to be introduced if we wish to have an NDAF. We may use this information 
to deduce a few interesting results on GRB models. 

NS-NS and BH-NS merger models, with (rout^^d) = (10,0.1) and (10,0.5) (see Popham et al. 

1999) , are well inside the NDAF zone and, according to our calculations, are capable of producing 
GRBs. However, this is only if the black hole is small (few Mq). If the black hole is larger than 

~ IOMq its Schwarszchild radius becomes too large and there is not enough "room" for an NDAF 
solution around it. Moreover, the neutron star in this case is swallowed whole by the BH and it 
is not tidally disrupted to create an accretion disk. On the other hand, as already noted, unless 
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the viscosity is much smaller than what we have assumed (which we consider unlikely), such disks 
cannot produce long bursts lasting hundreds or even tens of seconds. This suggests that NS-NS 
mergers and BH-NS mergers with smallish BH masses produce the class of short duration GRBs, 
but not the long duration GRBs. 

Other merger models, specifically the BH-WD and the BH-He star merger models, would 
appear not to be viable ORB engines. As the secondaries in these systems are not compact, they 
would form accretion flows with large values of rout- For instance, Popham et al. (1999) estimate 
^out ~ 3000 for a BH-WD binary and rout ~ 5000 for a BH-He star binary. At these radii, the 
accretion flow will be a very exteneded CDAF and hardly any mass will be accreted. Although 
the time scales of these models are consistent with long bursts, the extremely small value of rriacc 
suggests that these models do not produce GRBs of any kind. It is interesting to speculate what 
kind of observable events these binaries might produce (as undoubtedly they do merge in nature). 

All of the discussion so far is concerned with binary mergers, where we have imagined that a 
certain fixed amount of mass is instantaneously input into the accretion flow. The popular collapsar 
model (MacFadyen & Woosley 1999) corresponds to a different scenario in which mass is steadily fed 
over a period of time by fallback from the supernova explosion. MacFadyen and Woosley (1999) 
show that the time scale of the GRB is set by the physics of fallback rather than by accretion. 
Further, the time scales they obtain are consistent with observations of long GRBs. 

While the time scale may be set by fallback, the efficiency of the burst still depends on the 
nature of the post-fallback accretion. Let us ssume that fallback supplies mass at a certain injection 
rate rhi^j at a characteristic radius rout! the latter depends on the specific angular momentum of the 
material. Figure 3 shows the numerical results. As expected, efficient accretion, where most of the 
fallback material reaches the black hole, is possible only if rout is small and falls within the NDAF 
zone. If coUapsars have a distribution of rout ) then our calculations suggest that only those systems 
that have rout ~ lOOccI^^^mJ^ will make bursts. Systems with larger specific angular momentum, 
and hence larger rout; will form CDAFs and will eject most of the mass. Such systems may make 
very interesting supernova explosions, but if at all they make GRBs the bursts are likely to be very 
weak. 

We have assumed in this paper that the energy of a GRB is proportional to the total mass 
accreted on the central compact object. For NS-NS and NS-BH binaries, which give rise to short 
duration bursts, the disk mass is expected to be nearly constant and hence the total energy might 
be roughly the same for all bursts. However, for long bursts, which we assume result from the 
collapse of a massive star, the energy release depends on the angular momentum of the stellar core 
(which determines rout) and the mass of the stellar envelope that falls back on the collapsed core. 
Observations of long duration bursts indicate that the energy does not vary much from one burst 
to another (Panaitescu and Kumar 2001, Frail et al. 2001). This means that, for some reason, the 
total accreted mass on the central object depends only weakly on the properties of the progenitor 
star. 
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In our calculations we neglected photo-disintegration of nuclei, which Popham et al. (1999) 

and MacFadyen & Woosley (1999) show to be an important coolant of the accreting gas. If we 
include this effect, the boundary between CDAFs and NDAFs will move to somewhat higher values 
of Tout; perhaps by a factor of a few. However, it will not change our key conclusions. 

We are grateful to Eliot Quataert and Andrew MacFadyen for useful discussions, and Martin 
Rees for comments on the manuscript. We thank the organizers of the Jerusalem Winter School 
in Physics for hospitality while part of this research was carried out. This work was supported in 
part by NSF grant AST 9820686 and by a US-Israel BSF grant 98-00225. 
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A. Mass Accretion versus Ejection in CDAF Models 

We provide an alternate "derivation" of equation (10), and discuss the relevant uncertainties 
in the result. 

In a CDAF, there is a flow of energy from the inside to the outside as a result of convection. 
Let us write the convective luminosity as 

Lc = ecMinC^, 

where Min is the mass accretion rate onto the black hole. Since the convective luminosity is 
proportional to a (Quataert &; Gruzinov 2000; Narayan et al. 2000), we may write Cc = ar], with 
ri roughly a constant. Ball et al. (2001) computed a global model of a CDAF with a = 0.03, and 
found that Sc = 0.0045. This suggests that rj « 0.15. 

In our GRB model, there is practically no radiative emission from the gas. Therefore, the 
entire convective luminosity Lc must be converted into thermal and mechanical energy of the gas 
on the outside. This energy will cause the gas at r ~ rout to move to larger radii. 

How much energy do we need to evaporate gas from r = rout? Since the gas in a CDAF 
is hardly bound to the black hole — the Bernoulli parameter is close to zero or even positive 
(cf Quataert & Gruzinov 2000; Narayan et al. 2000) — very little energy is needed to drive a 
mass outflow to infinity. Let us write the energy carried away per unit mass of escaping gas as 
C,GM [rontRs- The rate at which mass is ejected is then given by 



LcToutRs 0.277a 

-(Wout = 

Thus we have 



,v ^crout-Ks 0.277a_i ■ 
= CGM = ^^^outMiu. 



.-1^-1 

ma Mout V 

Equation (10) of the main paper (derived by a different approach) has the same scaling, but with 
a coefficient equal to 14. 

By the Bernoulli parameter argument mentioned earlier, we expect C to be significantly less 
than unity; we would arbitrarily guess that ~ 0.1 — 0.3. Thus we expect bC,/r] ^ 10 (compared to 
14 in the main text). However, the coefficient is uncertain since C, could be much smaller than our 
assumed value, or conceivably much larger. In the latter case, the convective luminosity must be 
carried away by a small amount of gas that is accelerated to a speed much greater than the escape 
velocity. This possibility is not supported by the numerical simulations carried out to date. Our 
estimate of r] is also uncertain since it is based on a single global model calculated by Ball et al. 
(2001). A better estimate could be obtained from full-scale numerical simulations. 

The simulations reported by Igumenshchev & Abramowicz (2000) may be used to obtain a 
direct estimate of Min/Mout- For their model L (7 = 4/3, a = 0.03, rout ~ 7000), they estimate 
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Min/Mout ~ 0.003, whereas equation (10) of the present paper predicts a value of 0.007. The two 
estimates agree to within a factor of about 2, which is reassuring. On the other hand, for a = 0.1 
(Model I), they obtain Mjn/Mout ^ 0.02 which is much larger than our formula would predict. A 
possible explanation is that Model I does not strictly correspond to a CDAF. The flow appears 
to be a transition case where turbulent convection is replaced by a large scale circulation. In the 
case of MHD simulations (Stone Sz Pringle 2000) the mass outflow rate is again found to be much 
larger then the mass accretion rate on the central object. Therefore, the main results discused here 
should still be valid. 
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Fig. 1. — Contours of tacc in the rout (in units of Rs) ra^ (in solar masses) plane, for 7713 = 1 and 
a_i = 1. The contours are equally spaced in log. The lowest contour is at log(quantity) = —1.2, 
and succeeding contours are shifted by +0.2 up to a maximum value of +1. Two contours are 
highlighted: the contour at —1 is shown as a bold solid line and the contour at is shown as 
a bold dashed line. The three dotted lines correspond to three important transitions, computed 
with the analytical approximations of §§2,3. Starting from the right, the lines correspond to: (i) 
the transition from a CDAF to an NDAF (calculated via equation 20), (ii) the transition from a 
gas-pressure-dominated NDAF to a dcgcncracy-pressure-dominated NDAF (equation 31), and (iii) 
the transition from an NDAF that is optically thin to neutrino emission to one that is optically 
thick (equation 34). Our calculations are not reliable to the left of the leftmost dotted line since 
we have not allowed for optical depth effects in the model. 
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Fig. 2. — Contours of accretion efficiency ^ = macc/wzd in the Vout-n^d plane for = 1 and a_i = 1. 
The contour levels correspond to log^ = 0, —0.2, —0.4, —1 (solid lines) and —1.2, —1.8 (dotted 
lines). The white space on the left corresponds to the NDAF zone, where log^ = 0. 
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Fig. 3. — Contours of accretion efficiency ^ = m/rhinj in the Vout-'^nmi (in solar masses per second) 
plane. The contour levels and spacing are the same as in Fig. 2. 



